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IS.  ABSTRACT 


— — ^ An  investigation  has  been  conducted  to  determine  the  lift  hysteresis  loops 
on  a slender  elliptic  cylinder  undergoing  small  amplitude  angle  of  attack  oscil- 
lations about  some  mean  angle  of  attack.  The  results  obtained  in  this  investi- 
gation are  briefly  reviewed  in  this  final  report.  The  investigation  was  con- 
ducted in  three  phases.  In  the  first  phase  the  boundary  layer  over  the  cylinder 
was  assumed  to  be  entirely  laminar.  Simple  lift  hysteresis  loops  are  obtained 
in  this  case  and  these  hysteresis  loops  exhibit  some  of  the  same  general  charac- 
teristics as  the  lift  hysteresis  loops  for  airfoils  oscillating  in  pitch. 

The  second  phase  of  the  investigation  involved  a study  of  the  rature  un- 
steady boundary  layer  separation.  Analytical  solutions  were  obtained  for  the  un- 
steady laminar  boundary  layer  equations  for  several  flows  which  approach  separation. 
These  solutions  verify  the  Moore-Ratt-Sears  model  for  unsteady  boundary  layer 
separation,  at  least  for  the  case  of  upstream  moving  separation.  It  does  not 
appear  that  solutions  are  possible  for  the  case  of  downstream  moving  separation 
at  this  time. 

In  the  final  phase  of  this  investigation  an  attempt  was  made  to  determine 
the  lift  hysteresis  loops  for  the  case  where  the  boundary  layer  on  the  cylinder 
is  partially  turbulent.  It  was  found  that  a solution  to  this  problem,  by  the 
solution  technique  proposed,  is  not  possible  at  this  time. a 
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INTRODUCTION 

The  present  report  is  the  final  report  on  an  investigation  made  to  deter- 
mine the  unsteady  lift  hysteresis  on  a simple  body,  an  elliptic  cylinder,  os- 
cillating in  a uniform  stream.  The  purpose  of  the  investigation  was  to  pro- 
vide some  insight  into  the  problem  of  lift  hysteresis  or  dynamic  stall  which 
occurs  on  an  airfoil  oscillating  about  some  mean  angle  of  attack,  in  a uniform 
stream.  This,  in  turn,  is  a simplification  of  the  practical  problem  of  dynamic 
stall  encountered  by  a helicopter  blade  during  each  revolution  of  the  rotor. 

The  investigation  was  supported  by  the  U.  S.  Army  Research  Office  under 
grants  DA-ARO-D-31-124-72-G-134  and  DA  H C04-75-G-0007 , both  entitled  "An  In- 
vestigation of  Unsteady  Stall  Phenomenon  on  a Slender  Elliptic  Cylinder".  The 
work  on_  these  grants  was  carried  out  during  the  period  from  March  15,  1972 
through  September  30,  1977. 

This  final  report  provides  a brief  statement  of  the  problem  investigated, 
and  presents  the  results  obtained  and  conclusions  reached  during  each  phase  of 
the  investigation,  a list  of  reports  and  papers  which  resulted  from  the  inves- 
tigation and  a list  of  the  personnel  who  participated  in  the  investigation. 

STATEMENT  OF  THE  PROBLEM 

The  purpose  of  the  investigation  was  to  determine  the  variation  of  lift 
with  time  for  an  elliptic  cylinder  which  oscillates  with  an  angle  of  attack 
amplitude  Aa,  about  some  mean  angle  of  attack,  a^,  in  a uniform  stream  of 
velocity  U.  The  amplitude  of  the  oscillations,  Aa,  was  assumed  small,  in 
which  case  the  velocity  at  the  upper  edge  of  the  boundary  layer,  u^,  is  given 
by: 
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Here  R = /sin  n + 6 cos  n,  6 is  the  fineness  ratio  of  the  elliptic  cylinder, q 
is  a parametric  angle  related  to  the  distance  measured  along  the  cylinder  sur- 
face, x,  by  dx/dn  = -R/2,  l is  the  length  of  the  major  axis  of  the  ellipse,  oj 
is  the  frequency  of  oscillation  and  t is  time.  The  unsteady  circulation  is 
assumed  to  be  given  by 

T(t)  = Tq  + ri0  Aa  cos  wt  + r Aa  sin  wt  . 

The  values  of  T ^ and  must  be  determined  so  that  the  rate  at  which  vorticity 
is  shed  into  the  wake  is  equal  to  the  rate  of  change  of  the  circulation  about 
the  cylinder.  This  criterion  is  given  by: 

dT_l-2  1 2 , d^s  - dXs  ... 

dt  2 U6  " 2 -^5  + dt  " U6  dt  (2) 

where  the  upper  and  lower  bars  denote  conditions  at  separation  on  the  upper  and 
lower  surfaces  of  the  body,  respectively. 

The  problem  then  becomes  one  of  determining  the  values  of  T^q  and  which 
satisfy  equation  (2).  To  accomplish  this,  values  of  T^q  and  are  first 
assumed  and  the  boundary  layer  development,  over  both  the  top  and  bottom  surfaces 
of  the  elliptic  cylinder,  is  calculated  up  to  separation.  Once  the  locations 
of  the  upper  and  lower  separation  points  were  determined,  equation  2 can  be  used 
to  determine  if  the  proper  values  of  and  had  been  assumed.  If  the 
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assumed  values  of  T ^ and  T^  did  not  satisfy  equation  (2) , new  values  were  as- 
sumed and  the  process  was  repeated.  An  iteration  scheme  was  set  up  so  that  one 
could  rapidly  converge  on  the  proper  values  of  T and  T^. 

Initially  the  problem  was  formulated  for  the  case  where  the  boundary  layer 
on  the  elliptic  cylinder  was  laminar.  A later  phase  of  the  investigation  ex- 
plored the  case  where  the  boundary  layer  was  assumed  to  become  turbulent  at 
the  points  of  minimum  pressure. 

It  is  clear  from  equation  (2) , that  the  unsteady  lift  characteristics  de- 
pend upon  the  nature  of  the  moving  separation  points  on  the  body.  The  analysis 
discus  ove  it  was  assumed  that  the  separation  points  were  the  points  of 

van  on  the  upper  and  lower  surfaces  of  the  body.  It  is  well  known 

th  h. iiis  is  not  a correct  assumption.  For  the  case  of  small  amplitude,  low 

frequency  oscillations,  however,  it  is  probably  a reasonable  assumption. 
Nevertheless,  it  was  decided  to  investigate  the  nature  of  unsteady  separation 
as  part  of  this  project  in  an  effort  to  determine  the  proper  manner  in  which 
equation  (2)  should  be  evaluated.  An  effort  was  made  then,  as  part  of  this  in- 
vestigation, to  determine  the  nature  of  unsteady  laminar  boundary  layer  separation. 

RESULTS  AND  CONCLUSIONS 

Laminar  Boundary  Layer  Calculation  - Oscillating  Elliptic  Cylinder 

In  the  first  phase  of  the  investigation  the  calculations,  indicated  earlier, 
were  carried  out  for  the  case  where  the  boundary  layer  on  the  elliptic  cylinder 
was  assumed  to  be  entirely  laminar.  The  results  for  this  case  are  presented  in 
Reference  1. 

The  solution  technique  employed  in  the  laminar  boundary  layer  analysis  was 

L 


a momentum  integral  technique.  From  these  calculations  it  was  possible  to 
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determine  the  motion  of  the  upper  and  lower  separation  points  and  the  lift 
hysteresis  loops  for  the  oscillating  cylinder. 

The  calculated  lift  hysteresis  loop  direction  was  found  to  change  from 
counterclockwise  to  clockwise  as  the  mean  angle  of  attack  was  increased  from 
two  degrees  to  angles  of  attack  within  the  steady  state  stall  region.  In 
addition,  the  hysteresis  loops  become  larger  in  size  at  angles  of  attack 
corresponding  to  and  beyond  steady  state  stall.  Both  of  these  results  are  in 
general  agreement  with  results  obtained  experimentally  for  lift  hysteresis 
curves  of  airfoils  oscillating  in  pitch. 

The  same  calculations  were  carried  out  for  the  potential  flow  corres- 
ponding to  a stationary  elliptic  cylinder  with  a free  stream  which  oscillates 
in  angle  of  attack.  The  results  obtained  in  this  case  were  compared  with  these 
obtained  for  the  case  of  an  oscillating  cylinder.  These  calculations  show  that 
for  these  two  cases  the  lift  hysteresis  loops  were  similar  in  structure  but 
opposite  in  direction.  At  the  angle  of  attack  for  which  this  calculation  was 
carried  out,  the  hysteresis  loop  for  the  oscillating  body  in  a stationary 
stream  was  clockwise  while  the  hysteresis  loop  for  the  stationary  body  in  an 
oscillating  stream  was  counterclockwise.  The  potential  flows  for  the  two  sys- 
tems differ  only  by  a term  proportional  to  the  time  rate  of  change  of  the 
angle  of  attack;  a term  which  describes  the  influence  on  the  potential  flow 
of  the  added  mass  effect  normally  observed  when  a body  is  accelerated  in  un- 
steady flow.  Additional  results  obtained  in  this  phase  of  the  investigation  are 
presented  in  Reference  1. 

The  Nature  of  Unsteady  Boundary  Layer  Separation 

As  has  been  pointed  out  earlier,  the  criterion  which  determines  the  unsteady 
circulation  about  the  body,  and  thus  the  lift  on  the  body,  includes  both  the 
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velocity  at  the  boundary  layer  edge  at  separation  and  the  location  of  separa- 
tion (equation  2) . In  carrying  out  the  calculations  for  the  laminar  boundary 
layer  case,  described  above,  separation  was  assumed  to  occur  at  the  point  of 
vanishing  shear.  It  is  well  known  that  this  is  not  correct  but  the  error 
introduced  by  this  assumption  for  s lall  amplitude  oscillations  at  low  frequen- 
cies should  be  negligible.  Nevertheless,  it  was  decided  to  look  closely  at 
the  problem  of  unsteady  separation  to  determine  whether  or  not  it  is  possible 
to  determine  just  what  the  criterion  for  separation  in  unsteady  flow  is.  The 
results  of  this  phase  of  the  overall  investigation  are  presented  in  References  2, 

3 and  4. 

In  the  first  of  these  investigations  it  was  shown  that  for  a certain  class 
of  unsteady  flows,  namely  those  in  which  the  velocity  at  the  edge  of  the  bound- 
ary layer  is  a linear  function  of  x and  t,  the  flow  appears,  in  the  appropriate 
moving  coordinate  system,  as  a steady  flow  over  a wall  moving  at  constant 
velocity.  Solutions  to  the  problem  in  this  moving  coordinate  system  are  straight- 
forward using  well  developed  numerical  solution  techniques.  These  solutions, 
once  obtained,  can  be  transformed  back  into  the  fixed  coordinate  system,  where 
the  flow  is  unsteady,  and  used  to  infer  the  nature  of  unsteady  boundary  layer 
separation.  Solutions  were  obtained  for  the  case  where  the  boundary  layer  edge 
velocity,  in  the  moving  coordinate  system,  decreased  linearly  with  distance. 

The  solutions  obtained  tended  to  verify  the  Moore-Rott-Sears  model  for  unsteady 
separation  in  which  unsteady  separation  is  characterized  by  the  simultaneous 
vanishing  of  the  velocity  and  the  shear,  at  a point  with  the  boundary  layer, 
as  viewed  in  a coordinate  system  moving  with  separation.  Solutions  were  only 
obtained  for  the  case  corresponding  to  upstream  moving  separation. 

In  References  3 and  4 more  general  problems  were  considered.  It  was  shown 
that  the  method  of  semi-similar  solutions  was  well  suited  for  obtaining  solutions 
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to  the  unsteady  boundary  layer  equations  and  that  solutions  obtained,  using 
this  method,  for  flows  leading  to  separation  could  be  used  to  infer  the  nature 
of  unsteady  separation.  Reference  3 considered  a flow  which  is  an  unsteady 
variation  of  Howarth's  linearly  retarded  flow  while  Reference  4 considered 
an  unsteady  variation  of  the  Falkner  Skan  flows.  In  each  case  the  solutions 
verify  the  Moore-Rott-Sears  model  for  unsteady  boundary  layer  separation  in 
the  case  where  separation  is  moving  forward  along  the  body.  It  has  not  been 
possible  as  yet  to  obtain  solutions  for  the  case  where  separation  moves  back- 
ward along  the  body.  Williams  (Reference  5)  has  shown  that  in  the  cast  Jhere 
the  separation  point  is  moving  downstream  the  problem  is  still  parabolic,  but 
is  not  well  posed.  For  this  reason  it  has  not  been  possible,  as  yet,  to  ob- 
tain solutions  which  indicate  the  nature  of  downstream  moving  separation. 

The  method  of  semi-similar  solutions,  which  proved  quite  successful  in 
determining  the  nature  of  unsteady,  upstream  moving  separation,  was  extended 
to  the  problem  of  unsteady  turbulent  boundary  layers  in  Reference  6.  It  was 
shown  that  the  method  could  be  applied  to  certain  turbulent  boundary  layer 
flows.  The  question  of  unsteady  turbulent  boundary  layer  separation  was  not 
addressed,  however,  since  it  was  felt  that  the  algebraic  representation  for 
the  eddy  viscosity  employed  in  Reference  6 was  not  realistic  for  flows  which 
separate. 

Turbulent  Boundary  Layer  Calculations  - Oscillating  Elliptic  Cylinder 

In  the  third  and  final  phase  of  this  investigation  an  attempt  was  made  to 
determine  the  lift  hysteresis  loops  for  the  oscillating  elliptic  cylinder  in 
the  case  where  the  boundary  layer  on  the  cylinder  is  partially  turbulent.  The 
potential  flow,  given  by  equation  1,  may  be  written  simply  as 
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u.(x,t)  = u.  (x)  + u.  (x)  Aa  cos  u>t  + u_  (x)  Aa  sin  ait.  (3) 

6 <$0  62 

The  logical  method  of  solution  of  the  boundary  layer  equations  in  the  case 
where  the  external  velocity  is  of  this  form  is  to  expand  the  solution  in  a 
series  in  the  parameter  Aa.  The  first  term  in  such  a series  is  the  steady 
state  solution  for  the  elliptic  cylinder  at  the  mean  angle  of  attack,  a^. 

In  order  to  determine  the  vorticity  shed  into  the  wake  at  separation,  it  is 
necessary  to  solve  the  equations  for  each  term  in  the  series  up  to  the  separa- 
tion point.  In  the  laminar  case  this  is  no  problem  for  the  momentum  integral 
method  passes  through  separation  without  any  trouble.  In  the  turbulent  case, 
a number  of  different  solution  techniques  were  tried.  The  finite  difference 

technique  of  Cebeci  and  Smith and  the  integral  methods  of  Nash  and  Hicks^ 
(9) 

and  of  White  all  failed  to  predict  a definite  separation  point  for  the  steady 

state  case  (the  first  term  in  the  series).  Attempts  were  made  to  include  an 

interaction  effect  in  the  calculation  using  the  technique  devised  by  Dancey 

and  Pletcher^^,  but  this  technique  also  failed  to  predict  separation.  Finally, 

(8) 

the  method  of  Nash  and  Hicks v was  reformulated  using  the  velocity  profile 
given  by  Altstatt It  was  felt  that  this  velocity  profile  might  eliminate 
the  problems  encountered  because  it  does  not  contain  a logarithmic  singularity 
as  the  friction  velocity,  u , approaches  zero.  It  turns  out,  however,  that  this 

T 

method  also  runs  into  difficulty  as  separation  is  approached  and,  therefore, 
cannot  be  used  to  predict  separation. 

Apparently  the  problem,  in  the  case  of  turbulent  momentum  integral  methods, 
results  from  the  fact  that  the  product  of  the  boundary  layer  thickness,  6,  and 
the  friction  velocity,  u^,  appears  throughout  the  calculations.  As  separation 
is  approached  the  friction  velocity,  u^,  approaches  zero  but  the  boundary  layer 
thickness  becomes  large.  The  product,  therefore,  behaves  badly  unless  extreme 


care  is  used  in  the  calculation. 
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In  view  of  the  above  results  it  was  concluded  that  the  lift  hysteresis 
loops  could  not  be  obtained  at  this  time  for  the  case  where  the  boundary 
layer  is  turbulent.  This  problem  deserves  further  attention. 
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